ISSN 1364-0380 (on line) 1465-3060 (printed) 



1253 



Qeometry & Topology 
Volume 9 (2005) 1253-1293 
Published: 24 July 2005 




The colored Jones function is g— holonomic 

Stavros Garoufalidis 
Thang T Q Le 

School of Mathematics, Georgia Institute of Technology 
Atlanta, GA 30332-0160, USA 



Email: [stavro sOmath . gatech . eduj , [letuQinath . gatech . ed.u| 



URL: [http: /7www. math .gatech ■ edu/~stavros| 



Abstract 

A function of several variables is called holonomic if, roughly speaking, it is 
determined from finitely many of its values via finitely many linear recursion 
relations with polynomial coefficients. Zeilberger was the first to notice that 
the abstract notion of holonomicity can be applied to verify, in a systematic 
and computerized way, combinatorial identities among special functions. Using 
a general state sum definition of the colored Jones function of a link in 3~space, 
we prove from first principles that the colored Jones function is a multisum of a 
q-proper-hypergeometric function, and thus it is g-holonomic. We demonstrate 
our results by computer calculations. 
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1 Introduction 

1.1 Zeilberger meets Jones 

The colored Jones function of a framed knot /C in 3-space 

is a sequence of Laurent polynomials that essentially measures the Jones poly- 
nomial of a knot and its cables. This is a powerful but not well understood 
invariant of knots. As an example, the colored Jones function of the 0-framed 
right-hand trefoil is given by 



l/2-n/2 

J>c{n) = ^ — q-'^{l - - g^-") ... (1 



A:=0 



Here Jjcin) denotes the Jones polynomial of the 0-framed knot IC colored 
by the n-dimensional irreducible representation of 5(2, and normalized by 

Junknot(n) = _ ^-n/2)/(gl/2 _ g-l/2)_ 

Only a handful of knots have such a simple formula. However, as we shall see 
all knots have a multisum formula. Another way to look at the colored Jones 
function of the trefoil is via the following 3-term recursion formula: 

^n— 1 _j_ q4—4n g~™ (^^~^" q^~^^ ^—"in 

-^^("^ = gl/2(gn-l_g2-n) •^'^(^ ' ^) + q2-n _ gU-l ^^^^ " 2) 

with initial conditions: J/c(0) = 0, Jk{^) = 1- 

In this paper we prove that the colored Jones function of any knot satisfies a 
linear recursion relation, similar to the above one. For a few knots this was 
obtained by Gelca and his colleagues ^1 (In UH] a more complicated 

5-term recursion formula for the trefoil was established). 

Discrete functions that satisfy a nontrivial difference recursion relation are 
known by another name: they are g-holonomic. 

Holonomic functions were introduced by IN Bernstein [21 |3j and M Saito. The 
latter coined the term holonomic, that is a function which is entirely determined 
by the law of its differential equation, together with finitely many initial condi- 
tions. Bernstein used holonomic functions to prove a conjecture of Gelfand on 
the analytic continuation of operators. Holonomicity and the related notion of 
L'-modules are a tool in studying linear differential equations from the point 
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of view of algebra (differential Galois theory), algebraic geometry, and cate- 
gory theory. For an excellent introduction on holonomic functions and their 
properties, see |S] and |7j- 

Our approach to the colored Jones function owes greatly to Zeilberger's work. 
Zeilberger noticed that the abstract notion of holonomicity can be applied to 
verify, in a systematic and computerized way, combinatorial identities among 
special functions, and also [33l [28]. 

A starting point for Zeilberger, the so-called operator approach, is to replace 
functions by the recursion relations that they satisfy. This idea leads in a 
natural way to noncommutative algebras of operators that act on a function, 
together with left ideals of annihilating operators. 

To explain this idea concretely, consider the operators E and Q which act on 
a discrete function (that is, a function of a discrete variable n) /: N — > '^[q^] 
by: 

(Q/)(n) = <z"/(n) {Ef){n)=fin+1). 

It is easy to see that EQ = qQE , and that E, Q generate a noncommutative q- 
Weyl algebra generated by noncommutative polynomials in E and Q , modulo 
the relation EQ = qQE: 

A = Z[q^]{Q,E)/{EQ = qQE) 

Given a discrete function / as above, consider the recursion ideal If = {P G 
A \Pf = 0}. It is easy to see that it is a left ideal of the q-Weyl algebra. We 
say that / is q~holonomic iS If ^ 0. 

In this paper we prove that: 

Theorem 1 The colored Jones function of every knot is q -holonomic. 
Theorem n and its companion Theorem |2 are effective, as their proof reveals. 
Theorem 2 

(a) The E -order of the colored Jones function of a knot is hounded above hy 
an exponential function in the number of crossings. 

(b) For every knot /C there exist a natural number n(/C) , such that n(/C) initial 
values of the colored Jones function determine the colored Jones function of /C . 
In other words, the colored Jones function is determined by a Gnite list. n(/C) 
is bounded above by an exponential function in the number of crossings. 
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Computer calculations are given in Sectional In relation to (b) above, notice 
that the g-Weyl algebra is noetherian; thus every left ideal is finitely generated. 
The theorem states more, namely that the we can compute (via elimination) a 
basis for the recursion ideal of the colored Jones function of a knot. 

Let us end the introduction with some remarks. 

Remark 1.1 The colored Jones function can be defined for every simple Lie 
algebra q. Our proof of Theorem ^ generalizes and proves that the g -colored 
Jones function of a knot is g-holonomic (except for G2), see Theorem El below. 

Remark 1.2 The colored Jones function can be defined for colored links in 
3-space. Our proof of Theorem ^ proves that the colored Jones function of a 
link is g-holonomic in all variables, see Section f3. II 

Remark 1.3 It is well known that computing J{n) for any fixed n > 1 is a 
^P-complete problem. Theorem ^ claims that this sequence of #P-complete 
problems is no worse than any of its terms. 

Remark 1.4 The proof of Theorem ^ indicates that many statistical me- 
chanics models, with complicated partition functions that depend on several 
variables, are holonomic, provided that their local weights are holonomic. This 
observation may be of interest to statistical mechanics. 

1.2 Synonymous notions to holonomicity 

We have chosen to phrase the results of our paper mostly using the high-school 
language of linear recursion relations. We could have used synonymous terms 
such as linear g-difference equations, or g~holonomic functions, or D-modules, 
or maximally overdetermined systems of linear PDEs which is more common in 
the area of algebraic analysis, see for example [2l]. The geometric notion of D- 
modules gives rise to geometric invariants of knots, such as the characteristic 
variety introduced by the first author in The characteristic variety is 

determined by the colored Jones function of a knot and is conjectured to be 
isomorphic to the 5(2 (C) -character variety of a knot, viewed from the boundary 
torus. This, so-called AJ Conjecture, formulated by the first author is known 
to hold for all torus knots (due to Hikami, ^\), and infinitely many 2-bridge 
knots (due to the second-author, ^T]). 

Thus, there is nontrivial geometry encoded in the linear recursion relations of 
the colored Jones function of a knot. 
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1.3 Plan of the proof 

In Section [21 we discuss in detail the notion of a g-holonomic function. We give 
examples of (7~liolonomic functions (our building blocks), together with rules 
that create g-holonomic functions from known ones. 

In Section |31 we discuss the colored Jones function of a link in 3-space, using 
state sums associated to a planar projection of the link. The colored Jones 
function is built out of local building blocks (namely, 72 -matrices) associated 
to the crossings, which are assembled together in a way dictated by the planar 
projection. The main observation is that the i?-matrix is g-holonomic in all 
variables, and that the assembly preserves g-holonomicity. Theorem ^ follows. 
As a bonus, we present the colored Jones function as a multisum of a g-proper 
hypergeometric function. 

In Section 0] we show that the cyclotomic function of a knot (a reparametriza- 
tion of the colored Jones function, introduced by Habiro, with good integrality 
properties) is (7-holonomic, too. We achieve this by studying explicitly a change 
of basis for representations of 5l2- 

In Section 13 we give a theoretical review about complexity and computability of 
recursion relations of (7-holonomic functions, following Zeilberger. These ideas 
solve the problem of finding recursion relations of g-holonomic functions which 
are given by multisums of q proper hypergeometric functions. It is a fortunate 
coincidence (?) that the colored Jones function can be presented by such a 
multisum, thus we can compute its recursion relations. Theorem |21 follows. 

Section El is a computer implementation of the previous section, where we use 
Mathematica packages developed by A Riese. 

In Section 13 we discuss the g -colored Jones function of a knot, associated to a 
simple Lie algebra g. Our goal is to prove that the g-colored Jones function is 
g-holonomic in all variables (see Theorem El). In analogy with the g = s[2 case, 
we need to show that the local building block, the -matrix, is (/-holonomic in 
all variables. This is a trip to the world of quantum groups, which takes up the 
rest of the section, and ends with an appendix which computes (by brute-force) 
structure constants of quantized enveloping Lie algebras in the rank 2 case. 
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2 g— holonomic and g— hypergeometric functions 

Theorem ^ follows from the fact that the colored Jones function can be built 
from elementary blocks that are g-holonomic, and the operations that patch 
the blocks together to give the colored Jones function preserve q-holonomicity. 

IN Bernstein defined the notion of holonomic functions / : M^' — > C , [H El • 
For an excellent and complete account, see Bjork |^. Zeilberger's brilliant idea 
was to link the abstract notion of holonomicity to the concrete problem of algo- 
rithmically proving combinatorial identities among hypergeometric functions, 
see [23 OHj and also 28 . This opened an entirely new view on combinatorial 
identities. 

Sabbah extended Bernstein's approach to holonomic functions and defined the 
notion of a g-holonomic function, see [HJ and also [S]. 



2.1 g— holonomicity in many variables 

We briefly review here the definition of q -holonomicity. First of all, we need an 
r-dimensional version of the g-Weyl algebra. Consider the operators Ei and 
Qj for 1 < i, j < r which act on discrete functions /: N*" — > ^[q^] by: 

{Qif){ni,...,nr) = g""/(ni, . . . ,nr) 

{Eif){ni, . . . ,nr) = f{ni,...,ni-i,ni + l,ni+i,...,nr). 

It is easy to see that the following relations hold: 

QiQj — QjQi E{Ej — EjEj^ 

QiEj = EjQi fori / j EiQi = qQiEi 



(Relq 



We define the q~Weyl algebra Ar to be a noncommutative algebra with pre- 
sentation 

_ Z[q^^]{Qi,...,Qr,E^,...,Er) 

(Rel,) 
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Given a discrete function / with domain N*" or and target space a 7j[q^^- 
module, one can define the left ideal If in Ar by 

If:={PeAr\Pf = 0}. 

If we want to determine a function / by a finite list of initial conditions, it does 
not suffice to ensure that / satisfies one nontrivial recursion relation if r > 2. 
The key notion that we need instead is q-holonomicity. 

Intuitively, a discrete function /: N*' — > 'L\q^\ is g-holonomic if it satisfies 
a maximally overdetermined system of linear difference equations with polyno- 
mial coefficients. The exact definition of holonomicity is through homological 
dimension, as follows. 

Suppose M = Ar/I , where I is a left ^^-module. Let be the sub-space of 
Ar spanned by polynomials in Qi,Ei of total degree < m. Then the module 
Ar/I can be approximated by the sequence Fm/{Fm Ci I),m = 1,2,.... It 
turns out that, for m >> 1, the dimension (over the fractional field Q((7)) of 
Fm/ {Fm^I) is a polynomial in m whose degree d{M) is called the homological 
dimension of M . 

Bernstein's famous inequality (proved by Sabbah in the g-case, [S^) states 
that d{M) > r, if M 7^ and M has no monomial torsions, ie, any non-trivial 
element of M cannot be annihilated by a monomial in Qi,Ei. Note that the 
left Ar-module Mf := Ar ■ f — Ar/If does not have monomial torsion. 

Definition 2.1 We say that a discrete function / is g-holonomic if d{Mf) < r. 

Note that if d{Mj) < r, then by Bernstein's inequality, either Mf = or 
d{Mf) = r. The former can happen only if / = 0. 

Although we will not use in this paper, let us point out an alternative cohomo- 
logical definition of dimension for a finitely generated Ar module M . Let us 
define 

c(M) := min{i G N | Ext^^(M, A) / 0}. 

Then the homological dimension d{M) := 2r — c(M) equals to the dimension 
d{M) defined above. 

Closely related to Ar is the q-torus algebra %■ with presentation 
^ Z[q^']{Qf\...,Qt\Et\...,Et') 
(Rel,) 

Elements of % acts on the set of functions with domain Z*" , but not on the set 
of functions with domain N*". Note that %. is simple, but Ar is not. If / is a 
left ideal of %■ then the dimension of %■/! is equal to that of Ar/{I fl A^). 
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2.2 Assembling g— holonomic functions 

Despite the unwelcoming definition of g-holonomic functions, in this paper we 
will use not the definition itself, but rather the closure properties of the set of 
g-holonomic functions under some natural operations. 

Fact 

• Sums and products of g-holonomic functions are g-holonomic. 

• Specializations and extensions of g^-holonomic functions are g-holonomic. 
In other words, if /(ni, . . . , nm) is g-holonomic, the so are the functions 

g{n2, rim) ■=f{a, ^2, . . . , rim) 
and h{ni, . . . ,nm,nm+i) :=/(«!, . . ■ ,nm)- 

• Diagonals of q-holonomic functions are g-holonomic. In other words, if 
/(ni, . . . , rim) is g-holonomic, then so is the function 

g{n2, rim) ■= /(ra2, ^2, ns, . . . , n„). 

• Linear substitution. If /(ni, . . . , rim) is g-holonomic, then so is the func- 
tion, g{n'^, . . . ,n'j^i), where each n'j is a linear function of rij. 

• Multisums of g-holonomic functions are g-holonomic. In other words, if 
/(ni, . . . ,nm) is g-holonomic, the so are the functions g and h, defined 
by 

b 

g{a, b,n2,..., rim) ■= ^ /(^i, n2, . . . , rim) 

ni=a 
oo 

h{a, n2, . . . , rim) ■= ^ /(^i, n2, . . . , rim) 

ni=a 

(assuming that the latter sum is finite for each a). 

For a user- friendly explanation of these facts and for many examples, see j351l33| 
and I2H]- 

2.3 Examples of g— holonomic functions 

Here are a few examples of g-holonomic functions. In fact, we will encounter 
only sums, products, extensions, specializations, diagonals, and multisums of 
these functions. In what follows we usually extend the ground ring to 
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the fractional field Q{g^^^), where D is a positive integer. We also use v to 
denote a root of g, = q. 

For n,k £ Z, let 

r ^ n n 



{n}k 



n 




k 





i=l i=l 

'Ui=i{n-i + l}, iik>0 
if A; < 

if A; < 

The first four functions are g-holonomic in n, and the last two, as well as the 
delta function 6n,k, are (^-holonomic in both n and k. 

2.4 g— hypergeometric functions 

Definition 2.2 A discrete function /: W — > Qiq) is q-hypergeometric iff 
Eif/f eQ{q,q^\...,q^^) for alH = 1, . . . , r . 

In that case, we know generators for the annihilation ideal of /. Namely, let 
Eif/f = {Ri/Si)\Q.=q^i for Ri,Si € Z[q,Qi, . . . ,Qr]. Then, the annihilation 
ideal of / is generated by SiEi — Ri . 

All the functions in the previous subsections are g-hypergeometric. 

Unfortunately, g-hypergeometric functions are not always g-holonomic. For 
example, {n,k) — > l/[n^ + /c^]! is g -hyper geometric but not g-holonomic. 
However, g-proper-hypergeometric functions are g-holonomic. The latter were 
defined by Wilf-Zeilberger as follows, |3S1 Sec. 3.1]: 

Definition 2.3 A proper q-hypergeometric discrete function is one of the form 

F(n, k) = n.(^;g)a.rx+b..k+c. ^A(n,k)gk (-^^ 

where As,Bt G IK = Q{q) , as,ut are integers, bs,ks are vectors of r integers, 
A{n, k) is a quadratic form, CsjWs are variables and ^ is an r vector of elements 
in K. Here, as usual 

n — 1 

:=ll{l-Aq^). 

i=0 
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3 The colored Jones function for 
3.1 Proof of Theorem [T] for links 

We will formulate and prove an analog of Theorem ^ (see Theorem El below) 
for colored links. Our proof will use a state-sum definition of the colored Jones 
function, coming from a representation of the quantum group f/g(sl2), as was 
discovered by Reshetikhin and Turaev in |29| I32j . 

Suppose L is a framed, oriented link of p components. Then the colored Jones 
function Jl : W ^ = Z 

j^±i/2j defined using the representa- 

tions of braid groups coming from the quantum group Uq{5l2) ■ 

Theorem 3 The colored Jones function Jl is q-holonomic. 

Proof We will present the definition of Jl in the form most suitable for us. 
Let V{n) be the n-dimensional vector space over the field Q{v^^'^) with basis 
{eo, ei, . . . , e„_i}, with V{0) the zero vector space. 

Fix a positive integer m. A linear operator 



A: V{ni) (g) • • • V{nm) ^ V{n[) ® • • • ® F(n^) 
can be described by the collection 



where 



bi<n[,...,bm<n'^ 



We will call (oi, . . . , Om, 5i, . . . , coordinates of the matrix entry Aa\'^V.'.',am 

of A, with respect to the given basis. 

The building block of our construction is a pair of functions f±: Z^ ^ Z[t)^^/^] , 
given by 

f+{ni,n2;a,b,k) 



-(_l)V({ni-l-2a)(n2-l-26)+fc(fc-l))/2 b + k _ ^ ^ ^ _ 



f-{ni,n2;a,b,k) 

._ {{ni~l-2a-2k){n2-l-2b+2k)+k{k-l))/2 a + k , 
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The reader should not focus on the actual, cumbersome formulas. The main 
point is that: 

Fact 1 

• /+ and /_ are g -proper hyper geometric and thus g-holonomic in all 
variables. 

For each pair (ni,n2) G we define two operators 

B+{ni,n2),B-{ni,n2): V{ni)(^V{n2) V{n2)(^V{ni) 

by 

(^+(ni,n2))^'^^ := f+{ni,n2;a,b,c - b) 5c-b,a~d, 

(^_(ni,n2))^'^^ := U{ni,n2; a,b,b - c) 6c-b,a-d, 

where 6x,y is Kronecker's delta function. Although the coordinates (a, b, c, d) 
of the entry ;S±(ni, 712))^'^ of the operators B±{ni,n2) are defined for < 
a,b < rii and < c, d < n2 , the above formula makes sense for all non-negative 
integers a, b,c,d. This will be important for us. The following lemma is obvious. 

Lemma 3.1 The discrete functions B±{ni,n2)'^'l are q-holonomic with re- 
spect to the variables (ni, n2, a, b, c, d) . 

If we identify V{n) with the simple n~dimensional C/g(s[2)-module, with ei,i = 
0, . . . ,n — 1 being the standard basis, then B+{ni,n2),B-{ni,n2) are respec- 
tively the braiding operator and its inverse acting on V{ni) V{n2) ■ This fact 
follows from the formula of the i? -matrix, say, in jl7i Chapter 3]. In particular, 
B-{ni,n2) is the inverse of ;B+(ni, 712) . If one allows a, b, c, d in jB±(ni, 71,2)^'^ to 

run the set N, then ^±(ni, n2)|j'^ define the braid action on the Verma module 
corresponding to V{ni),V{n2) ■ 

Let Bm be the braid group on m strands, with standard generators cti, cXm-i : 

I 1 \/ t I 
1 i i+l 

For each braid /? € Bm and (ni,... ,nm) € N*", we will define an operator 
r(/3) = r(/3)(ni, . . . ,77^), 

r(/3) : V{ni) ® • • • ® V{nm) ^ • • • V^(n^(™)), 
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where (3 is the permutation of {!,... ,m} corresponding to (5. The operator 
r(/9) is uniquely determined by the fohowing properties: For an elementary 
braid crj, we have: 

T{af^) = id®*-i ^B±{n,, id®— . 

In addition, if /3 = f3'/3", then r(/3) := r(/5')r(/5") • It is well-known that r(/3) 
is well-defined. 

From Fact and Lemma l3. II it follows that 

Lemma 3.2 For any braid (i € Bm, the discrete function r(/3)(ni, . . . ,nm), 
considered as a function with variables ni , . . . , and all the coordinates of 
the matrix entry, is q-holonomic. 

Let K be the linear endomorphism of V{ni) ® ■ ■ ■ ® V{nm) defined by 

The inverse operator is well-defined. 

Corollary 3.3 For any braid P € Bm, the discrete function 

f{P) ■.= T{(3){nu...,nm)xK-^ 
is q~holonomic in ni, . . . , rim all all of the coordinates of the matrix entry. 

In general, the trace of f(/3) is called the quantum trace of t(/3). Although the 
target space and source space maybe different, let us define the quantum trace 
of r(/3)(ni, . . .,nm)) by 

tr,(/3)(ni,...,n^) := T(/3)(ni, . . . , n^)^;-^-. 

l<i<m 0<ai<ni 

It follows from Fact that trq(/3)(ni, . . . ,nm) is g-holonomic in ni, . . . ,nm- 
Restricting this function on the diagonal defined by rii = n^^^i = 1, . . . , m, we 
get a new function Jp of p variables, where p is the number of cycles of the 
permutation /3. 

Suppose a framed link L can be obtained by closing the braid (3. Then the 
colored Jones polynomial Jl is exactly Jg. Hence Theorem ^ follows. □ 

Remark 3.4 In general, Jicin) contains the fractional power g^/^. If K has 
framing 0, then Jic'{n) := JK{n)/[n] G See HOI. 
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Remark 3.5 There is a variant of the colored Jones function J^i of a colored 
link L' where one of the components is broken. If /? is a braid as above, let us 
define the broken quantum trace tr^ by 

tr;(/3)(ni,...,n„) := Yl ■ ■ ■ '^rn)a\::::aZ\ai=o- 

2<i<m 0<ai <ni 

Restricting this function on the diagonal defined by nj = uj^^ji = 1, . . . , m, we 
get a new function Jg/ of p variables, where p is the number of cycles of the 
permutation ^. 

If V denotes the broken link which is the closure of all but the first strand of 
j3 , then the colored Jones function Jj^' of L' satisfies J^' = J pi . 

If L denotes the closure of the broken link L', then we have: 

Jl = Jl' X [A] 
where A is the color of the broken component of L' . 

3.2 A multisum formula for the colored Jones function of a 
knot 

In this section we will give explicit multisum formulas for the 5(2 -colored Jones 
function of a knot. The calculation here is computerized in Section |SJ 

Consider a word w = a^^ . . . a^" of length m written in the standard generators 
0"!, . . . ,cTs_i of the braid group Bm with m strands, where = ±1 for all i. 

w gives rise to a braid /? E , and we assume that the closure of /? is a knot 
fC. Let IC' denote long knot which is the closure of all but the first strand of /?. 

A coloring of /C' is a tuple k = {ki, . . . , kc) of angle variables placed at the 
crossings of JC' . 

Lemma 3.6 There is a unique way to extend a coloring k of IC' to a coloring 
of the crossings and part-arcs of fC' such that 

• around each crossing the following consistency relations are satisfied: 

b+k a-k b-k a+k 




a b a b 

• The color of the lower-left incoming part-arc is . 
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Moreover, the labels of the part-arcs are linear forms on k . 

Proof Start walking along the long knot starting at the incoming part-arc. 
At the first crossing, whether over or under, the label of the outgoing part-arc 
is determined by the label of the ending part-arc and the angle variable of the 
crossing. Thus, we know the label of the outgoing part-arc of the first crossing. 
Keep going. Since fC' is topologically an interval, the result follows. □ 

For an example, see Figure E 




Figure 1: A word w, the corresponding braid /3, its long closure /C', and a coloring of 
/C' 

Fix a coloring of /C' determined by a vector k. Let 6i(k) for i = l,...,m 
denote the labels of the top part-arcs of p. Let Xj(k) and yj(k) denote the 
labeling of the left and right incoming part-arcs at the ith crossing of /C' for 
j = 1, . . . ,c. According to Lemma 6i(k),a;j(k) and yj(k) are linear forms 
on k. 

It is easy to see that 

k>0 

where 

m c 
i=2 j=l 

is a g-proper hypergeometric function. Remark 13.51 then implies that 

Proposition 3.7 The colored Jones function of a long knot K,' is a multisum 
of a q -proper hypergeometric function: 

JjC'{n) = ^F^(n,k). 

k>0 

Geometry & Topology, Volume 9 (2005) 



The colored Jones function is q-holonomic 



1267 



Remark 3.8 If a long knot /C' is presented by a planar projection D with c 
crossings (which is not necessarily the closure of a braid), then similar to the 
above there is a g-proper hypergeometric function F£,{n^ k) of c + 1 variables 
such that Jjc'in) = Ylkyo ^Din,\s.) . Of course, depends on the planar 
projection. Occasionally, some of the summation variables can be ignored. 
This is the case for the right hand-trefoil (where the multisum reduces to a 
single sum) and the figure eight (where it reduces to a double sum). 

D Bar-Natan has kindly provided us with a computerized version of Proposition 

EZl p. 



4 The cyclotomic function of a knot is g— holonomic 

Habiro ^S] proved that the colored Jones polynomial (of sl2) can be rearranged 
in the following convenient form, known as the cyclotomic expansion of the 
colored Jones polynomial: For every 0-framed knot /C, there exists a function 

oo 

such that JK{n) = Cjc{k)S{n, k), 

k=l 



where S{n^k) := {n + k — l}2k-i/ {v — v 



V 



-1 



Note that S{n,k) does not depend on the knot fC. Note that J is determined 
from C and vice-versa by an upper diagonal matrix, thus C takes values in 
Q{q). The difficult part of Habiro's result is Cfc takes values in The 
integrality of the cyclotomic function is a crucial ingredient in the study of 
integrality properties of 3-manifold invariants, |15| . 

Theorem 4 The cyclotomic function Cjc- N — > of every knot K, is 

q-holonomic. 

Proof Habiro showed that Cic{n) is the quantum invariant of the knot IC with 
color 



P"{n) :-- 



{2n - l}2n-2 
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where V{n) is the unique n-dimensional simple s[2-module, and (retaining 
Habiro's notation with a shift n — > n — 1) P"{n) is considered as an element 
of the ring of 5(2~modules over Q(t^) . 

Using induction one can easily prove that 

n 

P"{n)=y^R{n,k)V{k), 

where R{k, n) is given by 

R{n,h) = (-1 



k=l 



\n-k {2^} 



2n 
n — k 



{2n - l}![2n] 

We learned this formula from Habiro 15 and Masbaum '25^. Since 

k 

and -R(n, A;) is g-proper hypergeometric and thus g-holonomic in both variables 
n and fe, it follows that Cyc is g-holonomic. □ 



5 Complexity 

In this section we show that Theorem ^ is effective. In other words, we give a 
priori bounds and computations that appear in Theorem [21 

5.1 Finding a recursion relation for multisums 

Our starting point are multisums of g-proper hypergeometric functions. Recall 
the definition 12.31 of a g -proper hypergeometric function F(n,k) from Section 
im and let G denote 

G{n) :=^F(n,k) 

k>0 

throughout this section. 

With the notation of Equation Wilf-Zeilberger show that: 
Theorem 5 (|SS1 Sec.5.2]) 

(a) F(?7-,k) satisfies a k-free recurrence relation of order at most 

r! 
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where B = maxs^f{|bs[, |vf|, [as|, \ut\} + max^^ 

,v \0'^,u\ where CLf^^u sre the coefR- 

cients of the quadratic form A . 

(b) Moreover, G{n) satisfies an inhomogeneous recursion relation of order at 
most J*. 

Let us briefly comment on the proof of this theorem. A certificate is an operator 
of the form 

r 

P{E, Q) + Y,{E^ - l)Ri{E, Ei,...,Er,Q,Qi,...,Qr) 

i=l 

that annihilates F{n, k) , where P and Ri are operators with P a polynomial 
in E,Q, with P ^ 0. Here E is the shift operator on n, Ei (for i = 1, . . . , r) 
are shift operators in ki , and Q is the multiplication operator by q'' and Qi (for 
i = 1, . . . , r) are the multiplication operator by q^"- , where k = (ki, . . . ,kj.) . 

The important thing is that P{E, Q) is an operator that does not depend on 
the summation variables k. A certificate implies that for all (n, k) we have: 

r 

P(E',(5)F(n,k) + ^^{Gi{n,ki,... ,ki_i,ki + 1, fcj+i, . . . , /c^) - 
1=1 

Gi{n, ki, . . . , ki-i, ki, kij-i, . . . , kj.)) = 0, 

where Gi{n, k) = RiF{n, k) . Summing over k > 0, it follows that G{n) satisfies 
an inhomogeneous recursion relation PG = error(n) . Here error(n) is a sum of 
multisums of g-proper hyper geometric functions of one variables less. Iterating 
the process, we finally arrive at a homogeneous recursion relation for G. 

How can one find a certificate given F{n, k) ? Suppose that F satisfies a k-/ree 
recursion relation AF = 0, where A = A{E, Q, Ei, . . . , Ej.) is an operator that 
does not depend on the Qi . Then, evaluating A at Ei = . . . Er = I, we obtain 
that 

r 

A = A{E, Q, !,...,!)+ J2iE^ - l)R^{E, Q,Ei,...,Er) 

i=l 

is a certificate. 

How can we find a k-free recursion relation for F? Let us write 

where 5 is a finite set, j = {ji,...,jr), E^ = E^^ ...Ep', and crij{Q) are 
polynomial functions in Q with coefficients in Q(g); see [SHI. The condition 
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AF = is equivalent to the equation {AF)/F = 0. Since F is g-proper 
hypergeometric, the latter equation is the vanishing of a rational function in 
Qi, . . . ,Qr- By cleaning out denominators, this is equivalent to a system of 
linear equations (namely, the coefficients of monomials in Qi are zero), with 
unknowns the polynomial functions Cj j. For a careful discussion, see [30]. As 
long as there are more unknowns than equations, the system is guaranteed to 
have a solution. estimate the number of equations and unknowns in terms 
of -F(n, k), and prove Theorem |S1 

Wilf-Zeilberger programmed the above proof, see [23- As time passes the 
algorithms get faster and more refined. For the state-of-the-art algorithms and 
implementations, see |261 127j and [30], which we will use below. 

Alternative algorithms of noncommutative elimination, using noncommutative 
Grobner basis, have been developed by Chyzak and Salvy, ^8,. In order for have 
Grobner basis, one needs to use the following localization of the g-Weyl algebra 

„ _ Q{q,Qi,...,Qr){Eu...,Er) 
(Rel,) 

and Grobner basis jSI- 

In case r = 1 , ;Bi is a principal ideal domain [71 Chapter 2, Exercise 4.5]. In that 
case one can associate an operator in Bi (unique up to units) that generates 
that annihilating ideal of G{n) . For a conjectural relation between this operator 
for the 5(2 -colored Jones function of a knot and hyperbolic geometry, see 

Let us point out however that none of the above algorithms can find generators 
for the annihilating ideal of the multisum G{n) . In fact, it is an open problem 
how to find generators for the annihilating ideal of G{n) in terms of generators 
for the annihilating ideal of F{n, k) , in theory or in practice. We thank M 
Kashiwara for pointing this out to us. 



5.2 Upper bounds for initial conditions 

In another direction, one may ask the following question: if a g-holonomic 
function satisfies a nontrivial recursion relation, it follows that it is uniquely 
determined by a finite number of initial conditions. How many? This was 
answered by Yen, |33]. If G is a discrete function which satisfies a recursion 
relation of order J*, consider its principal symbol a{q, Q) , that is the coefficient 
of the leading E-teira.. The principal symbol lies in the commutative ring 
Z[g^,(5^] of Laurrent polynomials in two variables q and Q. For every n, 
consider the Laurrent polynomial <j{q,q"') € lL\q^\. If cr{q,q"') ^ for all n, 
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then G is determined by J* many initial values. Since cr{q,Q) ^ 0, it follows 
that cr{q,q"') ^ for large enough n. In fact, in 34, Prop. 3.1] Yen proves that 
a{q, q^) 7^ if n > degq{a) , then a{q, q^) ^ 0, where the degree of a Laurrent 
polynomial in q is the difference between the largest and smallest exponent. 
Thus, G is determined by J** := J* + degg((7) initial conditions. 

Yen further gives upper bounds for deg^ (cr) in terms of the g-hypergeometric 
summand, see j341 Thm.2.9] for single sums. An extension of Yen's work to 
multisums, gives a priori upper bounds J** in terms of the g-hypergeometric 
summand. These exponential bounds are of theoretical interest only, and in 
practice much smaller bounds are found by computer. 



5.3 Proof of Theorem El 



Theorem [^follows from Proposition 13 . 71 together with the discussion of Sections 
OandO □ 

Our luck with the colored Jones function is that we can identify it with a 
multisum of a g-proper hyper geometric function. Are we really lucky, or is 
there some deeper explanation? We believe that there is a underlying geomet- 
ric reason for coincidence, which in a sense explains the underlying geometry 
of topological quantum field theory. We will postpone to a later publication 
applications of this principle to Hyperbolic Geometry; 



6 In computer talk 



In this section we will show that Proposition 13.71 can be implemented by com- 
puter. 

For every knot, one can write down a multisum formula for the colored Jones 
function, where the summand is q-hypergeometric. Occasionally, this multisum 
formula can be written as a single sum. There are various programs that can 
compute the recursion relations and their orders for multisums. In maple, one 
may use qEKHAD developed by Zeilberger [2H1- In Mathematica, one may use 
the qZeil .m and qMultiSum.m packages of RISC developed by Paule and Riese 

[2niizainni. 
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6.1 Recursion relations for the cyclotomic function of twist 
knots 



The twist knots Kp for integer p are shown in Figure EJ Their planar projec- 
tions have 2\p\ + 2 crossings, 2\p\ of which come from the full twists, and 2 
come from the negative clasp. 



P 

full twists 




Figure 2: The twist knot Kp, for integers p. For p = — 1, it is the Figure 8, for p = 
it is the unknot, for p = 1 it is the left trefoil and for p = 2 it is the Stevedore's ribbon 
knot. 



Masbaum, j25j . following Habiro and Le gives the following formula for the 
cyclotomic function of a twist knot. Let c{p, •) denote the cyclotomic function 
of the twist knot Kp. Rearranging a bit Masbaum's formula '25' Eqn.(35)], we 
obtain that: 

C{p,n) = (-l)"+lg"("+3)/2 

\q;q)n+k+i{q;q)n-k 



'^(^_l^k^k(k+l)p+k{k-l)/2^^2k+l _ {Q;Q)n 



k=0 

The above sum has compact support for each n. Now, in computer talk, we 
have: 



Mathematica 4.2 for Sun Solaris 

Copyright 1988-2000 Wolfram Research, Inc. 

— Motif graphics initialized — 
In[l] :=<< qZeil.m 



q-Zeilberger Package by Axel Riese — ©RISC Linz — V 2.35 (04/29/03) 

For p = —1 (which corresponds to the Figure 8 knot) the program gives: 

In[2] := qZeil[q-(n(n + 3)/2) (-l)"(n + k + 1) q~(-k(k + 1)) (q-(2k + 1) 
- l)qfac[q, q, n]/(qfac[q, q, n + k + 1] qfac[q, q, n - k] ) 
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q~(k(k - l)/2), {k, 0, Infinity}, n, 1] 
Out [2]= SUM[n] == SUM[-1 + n] 

which means that c(— 1, n) = c{—l,n — 1) in accordance to the discussion after 
[251 Thm.5.1] which states c(— l,n) = 1 for all n. 

For p = 1 (which corresponds to the left hand trefoil) the program gives: 

In[3]:= qZeil[q-(n(n + 3)/2) (-l)'(n + k + 1) q-(k(k + l))(q-(2k + 1) 

- l)qfac[q, q, n]/(qfac[q, q, n + k + 1] qfac[q, q, n - k] ) 
q-(k(k - l)/2), {k, 0, Infinity}, n, 1] 

1 + n 

Out [3]= SUM[n] == -(q SUM[-1 + n] ) 

which means that c(l,n) = — g"+^c(l,n — 1) in accordance to the discussion 
after ^ Thm.5.1] which states c(l,n) = (-l)"g'^('^+3)/2 for all n. 

Similarly, for p = 2 (which corresponds to Stevedore's ribbon knot) the program 
gives: 

In[4]:= qZeil[q-(n(n + 3)/2) (-l)'(n + k + 1) q~(2k(k + l))(q-(2k + 1) 

- l)qfac[q, q, n]/(qfac[q, q, n + k + 1] qfac[q, q, n - k] ) 
q~(k(k - l)/2), {k, 0, Infinity}, n, 1] 

Out [4]:= No solution: Increase order by 1 

which proves that c(2, n) satisfies no first order recursion relation. It does 
satisfy a second order recursion relation, as we find by: 

In[5]:= qZeil[q-(n(n + 3)/2) (-l)'(n + k + 1) q~(2k(k + l))(q-(2k + 1) 

- 1) qfac[q, q, n]/(qfac[q, q, n + k + 1] qfac[q, q, n - k] ) 
q~(k(k - l)/2), {k, 0, Infinity}, n, 2] 

2 + 2 n -1 + n 

Out [4]= SUM[n] == -(q (1 - q ) SUM [-2 + n] ) - 

1 + n n 2 n 

> q (1+q-q+q ) SUM[-1 + n] 

Thus, the program computes not only a recursion relation, but also the order 
of a minimal one. Experimentally, it follows that c{p, n) satisfies a recursion 
relation of order for all p. Perhaps one can guess the form of a minimal 
order recursion relation for all twist knots. 

Actually, more is true. Namely, the formula for c{p, n) shows that it is a q- 
holonomic function in both variables (p, n) . Thus, we are guaranteed to find 
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recursion relations with respect to n and with respect to p. Usually, recursion 
relations with respect to p for fixed n are called skein theory for the nth colored 
Jones function, because the knot is changing, and the color is fixed. 

Thus, g-holonomicity implies skein relations (with respect to the number of 
twists) for the nth colored Jones polynomial of twist knots, for every fixed n. 

For computations of recursion relations of the cyclotomic function of twist knots, 
we refer the reader to |12j . 

6.2 Recursion relations for the colored Jones function of the 
figure 8 knot 

The Mathematica package qMultiSum.m can compute recursion relations for 
(^-multisums. Using this, we can compute equally easily the recursion relation 
for the colored Jones function. Due to the length of the output, we illustrate 
this by computing the recursion relation for the colored Jones function of the 
Figure 8 knot. Recall from Equation for p = —1 and from the fact that 
c(— l,n) = 1 that the colored Jones function of the figure 8 knot is given by: 

oo 

fc=o 

In computer talk, 

In[6]:= qZeil[q~(nk) qf ac [q~ (-n-l) ,q~ (-1) ,k] qf ac [q~ (-n+1) ,q,k] , 
{k,0,lnf iiiity},n,2] 

-1 - n n 2 n 

q (q + q ) (~q + q ) 

Out [6]= SUMCn] == 

n 

-1 + q 

-2 + n -1 + 2 n 

(1 - q ) (1 - q ) SUM [-2 + n] 

> + 

n -3 + 2 n 

(1 - q ) (1 - q ) 

-2 - 2 n -1 + n 2 -1 + n 

> (q (1 - q ) (1 + q ) 
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4 4n 3 + 11 l + 2n 3 + 2n l + 3n 

> (q + q - q - q - q - q ) SUM[-1 + n])/ 

n -3 + 2 n 

> ((1 - q ) (1 - q )) 

This is a second order inhomogeneous recursion relation for the colored Jones 
function. A third order homogeneous relation may be obtained by: 

In [7]:= MakeHomRec[7., SUM[n]] 

2 + n 3 n 
q (-q + q ) SUM [-3 + n] 

Out [7]= 

2 n 5 2 n 
(q + q ) (-q + q ) 

-2-n 2 n 8 4n 6 + n 7 + n 3 + 2n 

> (q (q-q)(q+q -2q +q -q + 

4 + 2n 5 + 2ii l + 3n 2 + 3n 

> q - q + q - 2 q ) SUM[-2 + n]) / 

n 5 2 n 

> ((q + q ) (q - q )) + 

-1-n n 4 4n 2+n 3+n l+2n 

> (q (-q+q)(q+q +q -2q -q + 

2 + 2n 3 + 2n l + 3n 2 + 3n 

> q - q - 2 q + q ) SUM[-1 + n]) / 

1 + n n 
2 n 2 n q (-1 + q ) SUM[n] 

> ((q + q ) (-q + q )) + == 

n 2 n 

(q + q ) (q - q ) 

Of course, we can clear denominators and write the above recursion relation 
using the g-Weyl algebra A. Let us end with a matching the theoretical bound 
for the recursion relation from Section [5] with the computer calculated bound 
from this section. Using Theorem [51 it follows that the summand satisfies a 
recursion relation of order J* = 1^ + 1^ = 2. This implies that the colored Jones 
function of the Figure 8 knot satisfies an inhomogeneous relation of degree 2 as 
was found above. The program also confirms that the colored Jones function 
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of the Figure 8 knot does not satisfy an inhomogeneous relation of order less 
than 2. 

7 The colored Jones function for a simple Lie algebra 

Fix a simple complex Lie algebra q of rank I. For every knot /C and every 
finite-dimensional g -module V ^ called the color of the knot, one can define 
the quantum invariant JjciV) G Zlq^^^"^^], where D is the determinant of the 
Cartan matrix of g. Simple g-modules are parametrized by the set of dominant 
weights, which can be identified, after we choose fixed fundamental weights, with 
N^. Hence Jjc can be considered as a function J/c- ^ 

Theorem 6 For every simple Lie algebra other than G2 , and a set of fixed 
fundamental weights, the colored Jones function J/c- — > is q- 

holonomic. 

Hence the colored Jones function will satisfy some recursion relations, which, 
together with values at a finitely many initial colors, totally determine the 
colored Jones function Jjc . 

Remark 7.1 The reason we exclude the G2 Lie algebra is technical. Namely, 
at present we cannot prove that the structure constants of the multiplication of 
the quantized enveloping algebra of G2 with respect to a standard PBW basis, 
are (7-holonomic; see Eemark IA.3I We believe however, that the theorem also 
holds for G2. 

The proof occupies the rest of this section. We will define J/c using representa- 
tion of the braid groups coming from the -matrix acting on Verma modules 
(instead of finite-dimensional modules). We then show that the 72 -matrix is 
q-holonomic. The theorem follows from that fact that products and traces of 
g-holonomic matrices are g-holonomic. 

7.1 Preliminaries 

Fix a Cartan subalgebra f) of g and a basis {ai,...,a^} of simple roots for 
the dual space f)*. Let be the R-vector space spanned by ai,. . . ,ai. The 
root lattice Y is the Z-lattice generated by {ai, . . . ,01}. Let X be the weight 
lattice that is spanned by the fundamental weights Ai, . . . , A^. Normalize the 
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invariant scalar product (•, •) on f)j^ so that (a, a) = 2 for every short root a. 
Let D be the determinant of the Cartan matrix, then {x,y) € ^^Z for G X. 

Let Si, i = 1, . . . , ^, be the reflection along the wall a^. The Weyl group W is 
generated by Si,i = 1 . . . ,i, with the braid relations together with = 1 . A 
word u) = Sjj . . . Si^ is reduced if considered as an element of W , can not be 
expressed by a shorter word. In this case the length l(w) of the element w (zW 
is r. The longest element ujq in has length t = (dim(g) — i)/2, the number 
of positive roots of q . 

7.1.1 The quantum group U 

The quantum group U = Uqis) associated to g is a Hopf algebra defined over 
Q{v), where v is the usual quantum parameter (see ^3 ESI)- Here our v is 
the same as v of Lusztig and is equal to q of Jantzen ^7] , while our q is 
v'^ . The standard generators of U are Ea, Fa, Ka for a G {ai, . . . , a^} . For a 
full set of relations, as well as a good introduction to quantum groups, see [T7] . 
Note that all the 's commute with each other. 

For an element 7 G y, 7 = kiUi + • • • + k^a^, let := iC^i . . . ■ 

There is a Y -grading on defined by \Ea\ and = 0. If 

X is homogeneous, then 

K^x = i>("'l^l)xi^T,. 

Let U'^ be the subalgebra of U generated by the Ea , by the , and 
by the Ka ■ It is known that the map 

[X,X ,X ) ^ XX X 

is an isomorphism of vector spaces. 

7.1.2 Verma modules and finite dimensional modules 

Let A G X be a weight. The Verma module M(A) is a Z//-module with under- 
lying vector space U~ and with the action of U that is uniquely determined by 
the following condition. Here r/ is the unit of the algebra U~ : 

Ea ■ r] = for all a 

Ka-r] = u^"'^)// for ah a 

Fa • X = FaX for all a G {qi, . . . , a^}, x G U~ 
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If (A|aj) < for all i = then M(A) is irreducible. On the other hand 

if (A|aj) > for all i = !,...,£ (ie, A is dominant), then M(A) has a unique 
proper maximal submodule, and the quotient L{X) of M(A) by the proper 
maximal submodule is a finite dimensional module (of type 1, see |17j). Every 
finite dimensional module of type 1 of Z// is a direct sum of several L(X). 

7.1.3 Quantum braid group action 

For each fundamental root a € {o-i, • • • , a^} there is an algebra automorphism 
Ta'. U ^ U , as described in ^1 Chapter 8] . These automorphisms satisfy the 
following relations, known as the hraid relations, or Coxeter moves. 

If (a,/3) = 0, then T^Tp = TpT^. 

If {a, (3) = -1, then T^TpTa = TpTaTp. 

If (a,/3) = -2, then TaTpTaTp = TpTaTpTa. 

If (a,/3) = -3, then TaTpTaTfiTo^Tp = T/sTaTpTaTpTa . 

Note that the Weyl group is generated by Sa with exactly the above relations, 
replacing by Sa , and the extra relations = 1 . 

Suppose w = Si-^^ . . . Si^ is a reduced word, one can define 

T-w := Tq..^ . . . Tq,.^ . 

Then is well-defined: If w, w' are two reduced words of the same element 
in W , then = T^' . This follows from the fact that any two reduced presen- 
tations of an element of W are related by a sequence of Coxeter moves. 

7.1.4 Ordering of the roots 

Suppose w = Si^Si^ ■ ■ ■ Si^ is a reduced word representing the longest element 
ujQ of the Weyl group. For r between 1 and t let 

7r(w) := SijSjj . . . Sj^_i(ajJ. 

Then the set {7j,i = is exactly the set of positive roots. We totally 

order the set of positive roots by 71 < 72 < • • • < 74 . This order depends on 
the reduced word w, and has the following convexity property: If Pi,P2 are 
two positive roots such that Pi + /32 is also a root, then (3i + [32 is between Pi 
and P2. In particular, the first and the last, 71 and 7t, are always fundamental 
roots. Conversely, any convex total ordering of the set of positive roots comes 
from a reduced word representing the longest element of W . 
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7.1.5 PBW basis for , , and 

Suppose w = Sij^ . . . Si^ is a reduced word representing tlie longest element of 
W . Let us define 

fr{w) = Ta^^ Ta^^ . . . Ta^^__^ J. 

Tlien Icrl = 7r = — |/r| • (We drop w if there is no confusion.) 

If 'fr is one of the fundamental roots, 7,- = a G {ai, . . . , a^}, then er{w) = Ea, 
fr{w) = Fa (and do not depend on w). 

For t > J > i > 1 let U~[j, i] be the vector space spanned by /"V^-i^ • • • j 
for all nj,nj-i, . . . ,ni G N and let U^[i,j] the vector space spanned by 
e^^'G^^i • ■ ■ 6j ^ ; for all nj,nj-i, . . . , G N. It is known that = h{^[t, 1] 
and U+ =U+[l,t]. 

For n = (ni, . . . , nt) G N*, j = (ji, . . . je) G and m = (mi, . . . , mt) G N* let 
us define f", KJ and e"^ by 

f^{w) := /r . . . Kj := K,,a, . . . e^{w) := e^^ . . . e^. 

Then as vector spaces over Q(t')W^, and hi have Poincare-Birkhoff-Witt 
(in short, PBW) basis 

{f" I n G N*}, {e"" I m G N*}, {f^RJe'" ] n, m G N*, j G Z^} 

respectively, associated with the reduced word w . 

In order to simplify notation, we define S := N* x x N* , and Xo- := f" . 
Thus, 

{x. I a G 5} (4) 
is a PBW basis of U with respect to the reduced word w . 

7.1.6 A commutation rule 

For x,y gU homogeneous let us define 

[x, y]q := xy — v^^^^'^'^^^yx. 

Note that, in general, is not proportional to [x,y]q. 

An important property of the PBW basis is the following commutation rule, see 
[TK| . If i < j then [fi,fj]q belongs to - l,i + 1] (which is if j = z + 1). 
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It follows that L{ [j, i] is an algebra. This allows us to sort algorithmically non- 
commutative monomials in the variables /j. Also two consecutive variables 
always g-commute: [fi,fi+i]q = 0. 

Similarly, ifi < j then [cj, ej]g belongs to U~^[i+1, j — 1] (which is if j = i+l). 
It follows that U~^[i,j] is an algebra, and two consecutive variables always q- 
commute, [ej,ej+i]q = 0. 

7.2 g— holonomicity of quantum groups 

Suppose A: U ^ U is a linear operator. Using the PBW basis of U (see 
Equation Q)), we can present A by a matrix: 

a' 

with A'^' G Q{v) . We call {a, a') the coordinates of the matrix entry A"^' . 

Definition 7.2 We say that A is q-holonomic if the matrix entry A'^' , con- 
sidered as a function of (a, a') is (7-holonomic with respect to all the variables. 

A priori this definition depends on the reduced word w. But we will soon see 
that if A is g-holonomic in one PBW basis, then it is so in any other PBW 
basis. 



7.2.1 holonomicity of transition matrix 

Suppose Xa{w') is another PBW basis associate to another reduced word w' 
representing the longest element of W . Then we have the transition matrix 
between the two bases, with entries in Q{v) . The next proposition checks that 
the entries of the transition matrix are g-holonomic, by a standard reduction 
to the rank 2 case. 

Proposition 7.3 Except for the case of G2, the matrix entry M^' is q- 
holonomic witli respect to aU its coordinates. 

Proof Since any two reduced presentations of an element of W are related 
by a sequence of Coxeter moves, it is enough to consider the case of a single 
Coxeter move. Since each Coxeter move involves only two fundamental roots 
and all Tq, 's are algebra isomorphisms, it is enough to considered the case of 
rank 2 Lie algebras. For all rank 2 Lie algebras (except G2) we present the 
proof in Appendix. □ 



Geometry & Topology, Volume 9 (2005) 



The colored Jones function is q-holonomic 



1281 



7.2.2 Structure constants 

Recall the PBW basis {xg- 1 o" G 5} of the algebra U . The multiplication in U 
is determined by the structure constants c{a,a',a") G Q{v) defined by: 



We will show the following: 

Theorem 7 The structure constant c{cr,a' ,a") is q-holonomic with respect 
to all its variables. 

Proof will be given in subsection 17.4.51 

7.2.3 Actions on Verma modules are g holonomic 

Each Verma module M(A) is naturally isomorphic to U~ , as a vector space, 
via the map u ^ u-r]. Using this isomorphism we identify a PBW basis of 
with a basis of M(A), also called a PBW basis. If u €U, then the action of u 
on M(A) in a PBW basis can be written by a matrix with entries in Q(v) . 
We call (n, n') E N* x N* the coordinates of the matrix entry. 

Proposition 7.4 For every r with 1 < r < t, the entries of the matrices 
e^, are q~holonomic with respect to k, X, and the coordinates of the entry. 

This Proposition follows immediately from Theorem [7| and Fact 0. 

7.3 Quantum knot invariants 
7.3.1 The quasi- i?— matrix 

Fix a reduced word w representing the longest element of W . For each r, 1 < 




r <t, let 



©r ■.= ^ckfr(^e^, 



fcGN 



where 
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Here v.y = t)('>'l'>')/2 , and 



k 



[kw- = n 



v:y-v^' 



-1 • 



The main thing to observe is that Ck is g-holonomic with respect to k. Note 
that although Qr is an infinite sum, for every weight X £ X , the action of 
on M(A) (8) M(A) is well-defined. This is because the action of is locally 
nilpotent, ie, for every x £ M(A), there is k such that • x = 0. 

The quasi- R -matrix is: 

G := GtGt_i...Gi. 

We will consider G as an operator from M(A) M(A) to itself. There is a 
natural basis for M(A) (?) M(A) coming from the PBW basis of AI(X). 

Proposition 7.5 The matrix of G acting on M{X) in a PBW basis is q- 
holonomic with respect to all the coordinates of the entry and X . 

Proof It's enough to prove the statement for each Q^. The result for G,. 
follows from the fact that the actions of e^,/^ on M(A), as well as Ck, are 
g-holonomic in k and so are all the coordinates of the matrix entries, by Propo- 
sition [731 • □ 



7.3.2 The matrix and the braiding 

As usual, let us define the weight on M(A) by declaring the weight oi Fa-rj to 
be A — ^ Uiji , where n = (ni, . . . , nt) ■ The space M(A) is the direct sum of 
its weight subspaces. 

Let V: M(A) M(A) M{X) M(A) be the linear operator defined by 

P(x y) = u^d^l'Is'Dx (g) y. 

Clearly T> is g^-holonomic; it's called the diagonal part of the -matrix, which 
is R := QV. 

The braiding is B := Ra , where a{x0y) = y®x. Combining the above results, 
we get the following: 

Theorem 8 The entry of the matrix of the braiding acting on M(A) is q — 
holonomic with respect to all the coordinates and X. 

Remark 7.6 Technically, in order to define the diagonal part V, one needs to 
extend the ground ring to include a D-th root of v, since (A, fj.) , with X, fi € X , 
is in general not an integer, but belonging to 4Z. 
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7.3.3 (7 holonomicity of quantum invariants of knots 



First let us recall the definition of quantum knot invariant. 

Using the braiding B: M(A) M(A) one can define a representation of the 
braid group r: (M(A))®'" by putting 

r(cJi) :=id®^-i®-B(g)id®'^-*-^ 

Let p denote the half-sum of positive roots. For an element x (z U and an 
Z// -module V , the quantum trace is defined as 

tVg{x,V) :=tv{xK.2p,V). 

Suppose a framed knot JC is obtained by closing a braid f3 G Bm- We would 
say that the colored Jones polynomial is the quantum trace of t(/3) . However, 
since M(A) is infinite-dimensional, the trace may not make sense. Instead, we 
will use a trick of breaking the knot. Let K,' denote the long knot which is the 
closure of all but the first strand of [5. 

Recall that t{(5) acts on (M(A))®'". Let 

r(/5)(A)S;:::iGZb^'/''] 

be the entries of the matrix r(/3)(A) . We will take partial trace by first putting 
ni = n']^ = and then take the sum over all n2 = Ug , . . . , n^ = . The 
following lemma shows that the sum is actually finite. 



Lemma 7.7 Suppose ni = 0. There are only a finite number of collections of 
(n2, n3, . . . , Um) £ N^~^ such that 

is not zero. 



Proof Let M'(A) be the maximal proper Z//-submodule of M(A). Then 
L(A) = M(A)/M'(A) is a finite dimensional vector space. In particular it has 
only a finite number of non-trivial weights. Hence, all except for a finite number 
of fn,n G N*, are in M'(A). 

We present the coefficients 13±{X) graphically as in Figure [51 

Note that if {B± )^i,n^^ is not equal to 0, then fma can be obtained from fm 
by action of an element in U , and similarly, f^^ can be obtained from by 
action of an element in U . Thus if we move upwards along a string of the braid, 
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Figure 3: and 



,m2 
n2 



the basis element at the top can always be obtained from the one at the bottom 
by an action of U . 

Because the closure of /? is a knot, by moving around the braid one can get any 
point from any particular point. Because the basis element fo is not in M'(A) , 
we conclude that if 

is not 0, with ni = 0, then all the basis vectors fna, • • • , fnm ^^t in M'(A) , 
and there are only a finite number of such collections. □ 

R ecall that 2p is the sum of all positive roots. Let us define 

n2,.--nmeN*,ni=0, 

Prom Q'— holonomicity of r(/3)(A) it follows that Jx;'(A) is g-holonomic. Jjc'i^) 
is a long knot invariant, and is related to the colored Jones polynomial Jic of 
the knot /C by 

^c(A) = J^'(A) xdim,(L(A)), 

where L{X) is the finite-dimensional simple -module of highest weight A, and 
dimq(L(A)) is its quantum dimension, and is given by the formula 

iX+p,a) _ -(A+p,o) 

<i-«(i(^» = n _„-,.■.) ■ 

Since dimg(L(A)) is g-holonomic in A, we see that JfcW is g-holonomic. This 
completes the proof of Theorem El □ 

Remark 7.8 The invariant Jjc of long knots is sometime more convenient. 
For example, J/cW might contain fractional power of q, but (if /C' has framing 
0,) J^'(A) is always in see (20]. Also the function J)c' can be extended 

to the whole weight lattice. 
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7.4 Proof of Theorem [7| 

7.4.1 is (7 holonomic 

We will need the linear maps ra,r^: — > , as defined in ^| Chapter 6]. 
Their restriction to U" is uniquely characterized by the properties: 

ra{xy) = raix) y + v'-"'^''^^ X ra{y) r^(xy) = x r^(y) + t;^"'!^!) r^(x) y (5) 

and for any two fundamental roots a,/3, (see ^1 Eqn.(6.15.4)]) and 

1 — 7;2" 

rUF^) = r'^{F^) = 5^^p-—^F2-\ (6) 
where := gee [T7I Eqn.(8.26.2)]. 

Lemma 7.9 For a fixed q G {oi, . . . , a^}, tiie matrix entries of the operators 
{fa)^i ir'a)^ '■ ^~ are q~holonomic with respect to k and the coordinates 

of the matrix entry. Similarly, (r„)'^', {r'^)^ : — > lA^ are q-holonomic. 

Proof We give a proof for : U~ IA~ . The other case is similar. 

There is a reduced word w' = Si-^ . . . Si^ representing the longest element 
of W such that = a. Then w = Si^ ■ ■ ■ Si^Sa is another reduced word 
representing ljq, where a := —LJo{a). 

For the PBW basis of U~ associated with w it's known that 7^ = a, and thus 
ft = Fa. According to ^1 8.26.5], for every x in the algebra U~[t — 1, 1] , one 
has 

raix) = 0. 

Using Equations © and © and induction, one can easily show that for every 
X G U-[t - 1, 1], 

^ 1 _ ,,2nt-2i+2 
i=l " 

This formula, applied to x = f^li^ ■ ■ ■ /"^ , proves the statement. □ 

7.4.2 f/q(s(2) is g holonomic 

Lemma 7.10 Theorem\^ holds true for g = 5(2. 
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Proof The PBW basis for U is F'^K^E'^, with m,n G N and j G Z. First of 
all we know that 



2 = 



m 




n 


i 




i 



F^-'biKo^;2i-n-m,i)E^-\ 



where 



b{Ka;a,i) := 



v„ - V, 



-1 



m 
i 



is the usual 



Here for any root 7, one defines = v^"''"'^/'^ , and 
quantum binomial coefficient calculated with v replaced by Va ■ 
Hence 

00 

{F'^K^E'^){F'^' K^' E"^') = ^F"'+"''-'a{m,k,n,m',k',n',i)E''+'''-\ 

1=0 

where 

a{m, k, n, m' , k' , n', i) 

_ 2k{i-m')+2k'{i~n) 



n 
i 



[i\\b{K;2i - n-m',i) K 



k+k' 



m 
i 

The value of the function a is in Z[u^^][A'^"'^] . Consider the coefficient of K^' 
in a; one gets a function of m, n, k, m', n', k', i, r with values in which is 

clearly g-holonomic with respect to all variables. The lemma follows. □ 



7.4.3 E^,F^: U ^ U are g holonomic in k 

Proposition 7.11 For a fixed fundamental root a G {qi, . . . , a^}, tfie opera- 
tors E^,F^: U ^ lA of left multiplication are q~holonomic with respect to k 
and all the coordinates of the matrix entry. Similarly, the right multiplication 
hy E^,F^ are q-holonomic with respect to k and all the coordinates of the 
matrix entry 

Proof (a) Left multiplication by and right multiplication by E^. 

Choose w as in the proof of Lemma 17.91 Then ft = Fa and et = Eq, , and 
an element of the PBW basis has the form f^^^xKj^ye^^ . It's clear that left 
multiplication by F^ and right multiplication by Ea are Q'-holonomic. 

(b) Left multiplication by E^- 
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Choose a reduced word w = Si^^ . . . Si^ representing the longest element loq that 
begins with a: = a. We have the corresponding PBW basis fi,ei,i = 
1, . . . ,t with fi = Fa and ei = ■ Thus a typical element of the PBW basis 
has the form 

xF^^KpE^^y, (7) 



where x = /"* 
has r^(x) 
that 



■eT'- By [m 8.26.6], since x e U-[t,2], one 



0. Using formula |171 6.17.1], one can easily prove by induction 



i—ik 



1=0 



k 
i 



\x)E, 



k—i 



Using this formula one can move the E^ past x in the expression ((Jj), (there 
appear Tq, and Ka), then one moves E^ past Fa using the 512 case. The last 
step is moving past Kp is easy, since 

EaKp = V-^'^'^'^KpEa. 

Using Lemmas 17.91 and I7.1U1 we see that each "moving step" is g-holonomic. 
Hence we get the result for the left multiplication by E^- 

(c) Right multiplication by F^. 

The proof is similar. We use the same basis (jZj) as in the case b). For y, by 
Lemma 8.26 of |17| . one has Vaiy) = 0. Hence using induction based on the 
formula (6.17.2) of ^Jj one can show that 



i(n—i) 



(^"^ 



n 
i 



Fr'Ka\r'a)\y). 



Using this formula, and the results for fLemma l7.9() and 5I2 (Lemma 17. 1U|) 
we can move Fa to the right. □ 



7.4.4 Ta is (7 holonomic 

Proposition 7.12 For a fixed fundamental root a G {ai, . . . the hraid 

operator Ta'. U and its inverse T"^ are q-holonomic. 

Proof By Proposition 17.31 we can use any PBW basis. 

Choose a reduced word w' = Sj^ . . . representing the longest element loq 
that begins with a : Ofj^ = a . Then w = Si^ ■ ■ ■ Si^ Sa is another reduced word 
representing luq, where a is the dual of q: a = —LJo{a). 
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We use fr to denote fr{w), and to denote fr{w'). The relation between the 
two PBW basis of w and w is as follows: For 1 < r < t — 1 , 

Besides, = = f[, et = = e[ . 

We will consider the matrix entry of Tq, : U U where the source space is 
equipped with the PBW corresponding to w, while the target space with the 
PBW basis corresponding to w' . 

From |171 Chapter 8], recall that: 
Hence 

Ta{F^) = To^iK) = i-irv-"'^"'~^^F^K^. 

For a basis element Xa = ft^ ■ ■ ■ fi^Kpe^^ . . . e™* , we have 
Ta{xa) = da{nt,mt)K-''^E^' X {flr^-^ 

. . . {f[rK,^p{e[r^ . . . {eT'-' X i^r, 

where da(,nt,mt) := (-i)^t+mt^nt(nt-i)-mt(rm-i) _ 

The left or right multiplication by is (/-holonomic with respect to n and 
all the coordinates. The left multiplication by F"* , as well as the right mul- 
tiplication my is g-holonomic with respect to nt and all coordinates, by 
Proposition 17. 1 11 One then can conclude that is g-holonomic. 

The proof for : U —^U is similar. One should use the PBW basis of w' for 
the source, and that of w for the target. □ 



7.4.5 Proof of Theorem 13 

It is clear that for each j G , the operator K j : U ^ U of left multiplication 
is g-holonomic. 

Fix a reduced word w representing the longest element of W . It suffices to 
show that for each 1 < r < t the operators e^, : U (left multiplication) 
are g-holonomic with respect to all variables, including k. 

This is true if = E^ and fr = Fa, where a is one of the fundamental 
roots, by Proposition 17.111 But any Cr or fr can be obtained from Ea and 
Fa by actions of product of various T^. 's. Hence from Proposition 17. 12l we get 
Theorem d □ 
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A Appendix: Proof of Proposition [7.31 for A2 and B2 

In this appendix we will prove Proposition 17.31 for the rank 2 Lie algebras A2 
and B2 ■ We will achieve this by a brute-force calculation. 

First, let us discuss some simplification, due to symmetry. The transition ma- 
trix of U leaves invariant each of . On the transition matrix is 
identity. Hence it's enough to consider the restriction of the transition matrix 
in hl~ and lA^ . Furthermore, the Cartan symmetry (the operator r of |17j ^ 
reduces the case of IA~^ to that of hl~ . 

A.l The case of A2 

There are two fundamental roots denoted by a and /3. The set of positive roots 
is {a, f3,a + l3} . The reduced representations of the longest element of the Weyl 
group are w = S1S2S1 and w' = S2S1S2, where si = Sa and 82 = 8/3. 

The total ordering (see Section [?.1.4j) of the set of positive roots corresponding 
to w and w' are, respectively: 

(71,72,73) = (a, a + /?,/?) 
(71', 72', 73') = {f3,a + P,a). 
Notice that 7j/ = 73-4 for i = 1, . . . , 3. 

The PBW basis of U~ (see Section 17.1. 5|) corresponding to w and w' are, 
respectively: 

{frf2f! I m,n,p € N}, {f^^f^^f^, \ m,n,p G N}, 

where 

(/3, /2, /i) = {Fp,Ta{Fp) = -v[Fp, F^]g = FpF^ - vF^Fp, F^) 
(/3',/2',/i') = {Fa,Tp{Fa) = F^Fp-vFf3Fa,Fp). 
From explicit formulas of [231 section 5] it follows that: 

Lemma A.l The structure constants ofU'' , in the basis of w , is q~holonomic. 

Let us define a scalar product (•, •) on such that the PBW basis of w is an 
orthonormal basis. Since 

fm! fn' fV' \ ^ / fm' rn' fP' rm fn fn\ rm rn fP 

J3' J2' Jv — U3' J2' Jl' y Js J2 Jl )J3 J2 Jl 

m,n,p 
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Proposition 17.31 is equivalent to showing that 

r fm' fn' fp' fm m fn\ 
Us' /2' Tr ' J 3 J2 Jl I 

is g-holonomic in all variables m,n,p,m' ,n' ,p' . 

Since multiplication is g'-holonomic in the PBW basis of w (see Lemma lA.l|) . 
it suffices to show that 

/ fk fm fn fn\ 
\Ji' 5 Js J 2 Jl ) 

is g-holonomic in k, m, n,p for each i = 1, 2, 3. This is clear for i = 1 or z = 3, 
since fi/ = fs and fy = fi- As for f2' , an easy induction shows that 

oo 
k=0 

and the statement also holds true for i = 2. This proves Proposition 17.31 for 

A2. 



fk f-n—k fk 
JS J 2 J I ■ 



A. 2 The case of B2 

There are two fundamental roots denoted here by a and /3, where a is the 
short root. The set of positive roots is {a, /3,2a + (3,a + fi}. The reduced 
representations of the longest element of the Weyl group are w = S1S2S1S2 and 
w' = S2S1S2S1, where si = Sa and S2 = S/3- 

The total ordering of the set of positive roots corresponding to w and w' are, 
respectively: 

(71,72,73,74) = (a, 2a + /3, a + /?, /3) 
(71', 72', 73', 74') = (/3,a + /3,2a + /3,a). 
Notice that 7j/ = 74. j for z = 1, . . . , 4. 

The PBW basis of (see Section 17.1. 5p corresponding to w and w' are, 
respectively: 

{fifS'f2f?\ l,m,n,peN}, {fhf^f^>ff,\l,m,n,p&n}, 

where 

(/4,/3,/2,/i) = {Fp,FpFa-v^FaFp,^-vFaFpFa + ^^^^l^,Fa) 

if A' , /3' , /2' Jl') = {Fa , ^^-^ - vFaFpFa + ^ , FaFfS - FpFa ,Fp). 

It follows from [231 that: 
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Lemma A. 2 The structure constants ofU , in the basis of w , is q~holonomic. 

Let us define a scalar product (•, •) on U~ such that the PBW basis of w is an 
orthonormal basis. Then Proposition 17.31 is equivalent to 



is g-holonomic in all variables l,m,n,p,l'm' ,n' ,p' . 

Since multiplication is g-holonomic in the PBW basis of w (see Lemma lA.2|) . 
it suffices to show that 



is (7-holonomic in k,l,m,n,p for each i = 1,2,3,4. This is clear for z' = 1 or 
i' = 4 , since = fi and f^/ = /i . As for i' = 2 and i' = 3, the formula of 
[m Section 37.1] shows that 



and since Fa = f^' and = , the cases of i' = 2' and i' = 3' reduce to the 
cases of i' = 1' and i' = A' . This proves Proposition 17.31 for B2 ■ 

Remark A. 3 If Lemma lA . 1 1 holds for G2, then we can prove Proposition 17.31 
for G2. 



References 

[1] D Bar-Natan, Colored Jones. nb, Mathematica program, part of "KnotAtlas", 
April 2003. 

[2] IN Bernstein, Modules over a ring of differential operators. An investigation 
of the fundamental solutions of equations with constant coefficients^ Funkcional. 
Anal, i Prilozen. 5 (1971) 1-16, English translation: 89-101 MR0290097 

[3] I N Bernstein, Analytic continuation of generalized functions with respect to 
a parameter, Funkcional. Anal, i Prilozen. 6 (1972) 26-40, English translation: 
273-285 MR0320735I 

[4] J-E Bjork, Rings of differential operators, volume 21 of North-Holland Mathe- 
matical Library, North- Holland Publishing Co., Amsterdam (1979) MRMJIM, 

Geometry & Topology, Volume 9 (2005) 



a I' fm' fTi' fv' fl fm fn fn\ 
4'/3' J 2' Jy 5 /4/3 J 2 J I ) 



ak fl fm fn fP\ 
i' 5 74/3 T2 J I ) 





1292 



Stavros Garoufalidis and Thang TQ Le 



[5] A Borel, P-P Grivel, B Kaup, A Haefliger, B Malgrange, F Ehlers, 

Algebraic D -modules, Perspectives in Mathematics 2, Academic Press, Boston, 
MA (1987) MR882000 

[6] P Cartier, Demonstration "automatique" d'identites et fonctions hy- 
pergeometriques (d'apres D. Zeilberger), Asterisque (1992) Exp. No. 746, 3, 
41-91 MR1206064 

[7] S C Coutinho, A primer of algebraic D -modules, volume 33 of London Math- 
ematical Society Student Texts, Cambridge University Press, Cambridge (1995) 
MR1356713 

[8] F Chyzak, B Salvy, Non- commutative elimination in Ore algebras proves 
multivariate identities, J. Symbolic Comput. 26 (1998) 187-227 MRl 635242 

[9] S Garoufalidis, TTQ Le, D Zeilberger, The quantum MacMahon Master 
Theorem, [arXiv : math . qA/0303319| to appear in Proc. Nat. Acad. Sci. 

[10] S Garoufalidis, Difference and differential equations for the colored Jones func- 
tion, arXiv:math.GT/0306229 

[11] S Garoufalidis, On the characteristic and deformation varieties of a 
knot, from: "Proceedings of the CassonFest (Arkansas and Texas 2003)", 
Geom. Topol. Monogr. 7 (2004) 291-309 

[12] S Garoufalidis, X Sun, The C -polynomial of a knot, preprint (2005) 
larXiv : math . GT/0504305 

[13] R Gelca, Non-commutative trigonometry and the A-polynomial of the trefoil 
knot. Math. Proc. Cambridge Philos. Soc. 133 (2002) 311-323 MR1912404 

[14] R Gelca, J Sain, The noncommutative A-ideal of a (2, 2p + \) -torus knot 
determines its Jones polynomial, J. Knot Theory Ramifications 12 (2003) 187- 
201 IMR1967240i 

[15] K Habiro, On the quantum SI2 invariants of knots and integral homol- 
ogy spheres, from: "Invariants of knots and 3-manifolds (Kyoto, 2001)", 
Geom. Topol. Monogr. 4 (2002) 55-68 MR2002603 

[16] K Habiro, TTQ Le, in preparation 

[17] J Jantzen, Lectures on quantum groups, volume 6 of Graduate Studies in Math- 
ematics, American Mathematical Society, Providence, RI (1996) MR1359532 

[18] LI Korogodski, YS Soibelman, Algebras of functions on quantum groups. 
Part I, volume 56 of Mathematical Surveys and Monographs, American Mathe- 
matical Society, Providence, RI (1998) MR1614943 

[19] K Hikami, Difference equation of the colored Jones polynomial for torus knot. 
Internal. J. Math. 15 (2004) 959-965 MR2106155 

[20] TTQ Le, Jntegrality and symmetry of quantum link invariants, Duke Math. J. 
102 (2000) 273-306 <MRr74943ia 

[21] TTQ Le, The Colored Jones Polynomial and the A-Polynomial of Knots, 
..arXiv : math . GT/0407521. 



Geometry & Topology, Volume 9 (2005) 



The colored Jones function is q-holonomic 



1293 



[22] G Lusztig, Introduction to quantum groups, volume 110 of Progress in Mathe- 
matics, Birkhauser Boston Inc., Boston, MA (1993) MR1227098 

[23] G Lusztig, Quantum groups at roots of 1, Geom. Dedicata 35 (1990) 89-113 
MR1066560 

[24] B Malgrange, Equations differentielles a coefficients polynomiaux, volume 96 
of Progress in Mathematics, Birkhauser Boston Inc., Boston, MA (1991) 
MR1117227 

[25] G Masbaum, Skein-theoretical derivation of some formulas of Hahiro, 



Algebr. Geom. Topol. 3 (2003) 537-556 |MR1997328 



[26] P Paule, A Riese, A Mathematica q-analogue of Zeilberger's algorithm based 
on an algebraically motivated approach to q-hypergeometric telescoping, from: 
"Special functions, g-series and related topics (Toronto, ON, 1995)" , Fields Inst. 
Commun. 14, Amer. Math. Soc, Providence, RI (1997) 179-210 MR144868II 

[27] P Paule, A Riese, Mathematica software, available at: 

http : //www.risc .uni-linz . ac . at /re search/ comb inat /rise/ sof tware/qZeil/l 

[28] M Petkovsek, HS Wilf, Doron Zeilberger, A = B, A K Peters Ltd., 
Wellesley, MA (1996) MR1379802 

[29] N Yu Reshetikhin, V G Turaev, Ribbon graphs and their invariants derived 
from quantum groups. Comm. Math. Phys. 127 (1990) 1-26 MR1036112 

[30] A Riese, qMultisum-A package for proving q -hypergeometric multiple summa- 
tion identities, preprint (2002) 

[31] C Sabbah, Systemes holonomes d'equations aux q- differences, from: "Z?- 
modules and microlocal geometry (Lisbon, 1990)", de Gruyter, Berlin (1993) 
125-147 MR1206016 

[32] V G Turaev, The Yang-Baxter equation and invariants of links, Invent. Math. 
92 (1988) 527-553 MR939474 

[33] HS Wilf, D Zeilberger, An algorithmic proof theory for hypergeometric (or- 
dinary and "q ") multisum/ integral identities, Invent. Math. 108 (1992) 575-633 
MR1163239 

[34] L Yen, A two-line algorithm for proving q-hypergeometric identities, J. Math. 
Anal. Appl. 213 (1997) 1-14 MR1469359 

[35] D Zeilberger, A holonomic systems approach to special functions identities, J. 
Comput. Appl. Math. 32 (1990) 321-368 MR1090884 



Qeometry & Topology, Volume 9 (2005) 



